Possible uniform-flux chiral spin liquid states in the SU(3) ring-exchange model on the triangular 

lattice 
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We consider a SU(3) model with antiferromagnetic three-site ring exchanges, in addition to two-site ex- 
changes, on the triangular lattice. We first present numerical site-factorized state studies on the magnetic ordered 
states, which shows two different three-sublattice-ordered states, the antiferro-quadrupolar phase and the stan- 
dard 120° anti-ferromagnetic phase, along the axis of the strength of the three-site ring exchanges. We further 
study the model using slave-fermion mean field approaches in which we rewrite the exchange operators in terms 
of three flavors of fermions. At the mean-field level, we find the main competing trial states are the trimer state 
(triangular plaquette state), gapped uniform 7r/3-flux chiral spin liquid, and gapped uniform 27r/3-flux chiral 
spin liquid. The filled band of the 7r/3-fiux chiral spin liquid has Chern number +1, and that of the 27r/3-fiux 
chiral spin liquid has Chern number +2. We also give the effective Chern-Simons theory for each chiral spin 
liquid at the mean-field level. 



I. INTRODUCTION 

One of the goals in the studies of strongly interacting cold 
atom gases is that they can be used to simulate strongly cor- 
related systems^. The model systems can be engineered with 
a high degree of control and experimentally studied to reveal 
the profound nature of the phases, among which the Quantum 
spin liquid states (QSL)2^ are perhaps one of the most intrigu- 
ing phases. With cold atom gases, the iV-flavor fermionic 
Hubbard model on different lattices can be also realized. A 
model Hamiltonian to describe such systems is the A^-flavor 
fermionic Hubbard model^— 
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where a, j3 run over the different flavors, (jk) runs over pairs 
of nearest neighbors on the lattice, and j runs over all lattice 
sites. 

If we focus on the 1 /N filling, the system with generic A^ 
flavor of fermions can undergo metal-to-Mott insulator phase 
transition for sufficiently large repulsion U. The transition to 
a Mott insulator has been recently observed in (N = 2) spin- 
1/2 Hubbard models In this case, it is generally accepted 
that the ground state is well-captured by the usual antiferro- 
magnetic (AFM) Heisenberg model. On the other hand, when 
U is not very large compared with the hopping strength t, the 
ground state is not well-understood and it is possible that the 
strong charge fluctuations near the transition play an impor- 
tant role for stabilizing "gapless" spin liquid phases^— 

For N > 2ji2r— in the large U limit we also obtain 
Heisenberg-like (two-site exchange) Hamiltonian 



where Pjk is the two-site exchange operator, which per- 
mutes the fermions between two nearest-neighbor sites as 
Pjk\&, (3) = a), where the a, (3 represent the spin states 
at sites j and k. For N = 3, there has been numerical ev- 
idence on such a SU(3) Heisenberg model on the triangular 



lattice suggesting three-sublattice-ordered ground state in this 
regime^ 

However, at t/U ~ 0(1) in the SU(3) case, the two-site 
exchange Hamiltonian is not sufficient to capture the essential 
physics and higher ordered contributions should be consid- 
ered. If we perform perturbation studies on the SU(3) Hub- 
bard model at 1/3-filling, we will in fact obtain "ferromag- 
netic" (FM) three-site ring exchanges and four-site ring ex- 
changes and etc., if higher-order terms are included. In the 
FM three-site ring exchanges regime with four-site ring ex- 
changes included, in Ref. |22j we explored the magnetic or- 
dered phases using the site-factorized ansatz, which found 
three-sublattice ordered states and FM state in a large parame- 
ter regime; we also explored the non-magnetic ordered phases 
using the mean-field slave-fermion trial states. By qualita- 
tively arguing that the quantum fluctuations in linear flavor 
wave theory approximations destroy the three-sublattice or- 
dered states in strong ring exchanges, We conjectured that 
several possible gapless QSL were present due to the geomet- 
rical frustration along with the strong ring exchanges near the 
metal-Mott phase transition. 

For the AFM side of three-site ring exchange regime, at the 
first sight it seems unnatural to consider this regime. However, 
from the perspective of cold atom systems, recently the cold 
atom experiment by Struck et alM demonstrated a method to 
be able to add an artificial tunable gauge potential to the sys- 
tem. With the tunable gauge potential, it is possible to tune 
the sign of the three-site r ing exchanges from FM to AFM. 
In addition, recently, Ref. [24- did variational studies on the 
SU(3) model with three-site ring exchange (no higher-order 
terms such as four-site ring exchanges and etc.), and found, 
on the AFM side of the three-site ring exchanges, an inter- 
esting d x + idy chiral spin liquid state (CSL) with a gapless 
parton Fermi surface. However, the uniform-flux chiral spin 
liquid states^ - — QSL that break both parity and time-reversal 
symmetries, were not considered in the previous studies. In 
this paper, we focus on the SU(3) model with "AFM" three- 
site ring-exchange on the triangular lattice and we find, at the 
mean field level, the energies of uniform-flux CSL are lower 
than the interesting d x + id y CSL and suggest that they are 
stabilized in this model, at least at the mean field. 
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In order to obtain all possible ground states including mag- 
netic ordered and non-magnetic phases, in this work, we first 
study the magnetic ordered phases using the site-factorized 
ansatz3 We find that the phase diagram contains two differ- 
ent three-sublattice-ordered phases, the same to those found in 
Ref. I2II We further study this model using mean-field slave- 
fermion trial states focusing on the non-magnetic phases. Af- 
ter performing numerical full optimization of the trial energy, 
we consider three ansatz states-the uniform 7r/3-flux chiral 
spin liquid state CSL), the uniform 27r/3-flux chiral spin 
liquid state ($2^ CSL), and the trimer (plaquette) state. At the 
mean-field level, each gapped uniform-flux state breaks Time 
Reversal (TR) and possesses a finite Chern number, C = +1 

for the $ - CSL and C = +2 for the $ 2, CSL. 

3 — 

The paper is organized as follows. In Sec. |II]we define ex- 
plicitly the model Hamiltonian we will study. In Sec. Ill Al we 
use the site-factorized ansatz to study the magnetic ordered 
states. In Sec. lIIBI we use the slave-fermion representation to 
rephrase the SU(3) Hamiltonian in terms of the three-flavor 
fermionic Hamiltonian and perform the fermionic mean-field 
treatment of the model. In Sec. Ill CI we give the effective 
Chern-Simons theories for each CSL stabilized in this SU(3) 
ring exchange model at the mean-field level. In Sec. HTT1 we 
conclude with some discussions. 




FIG. 1 . The illustration of the triangular lattice. For each triangle, we 
label each site counterclockwise from 1 to 3. For each site, there are 
two triangles, a green backslashed up triangle and a blue forward- 
slashed down triangle, associated with it. 
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II. SU(3) MODEL WITH RING EXCHANGE TERMS 

The model Hamiltonian we consider is 

•* A 



FIG. 2. The phase diagram using the site-factorized states. There 
are two different three-sublattice-ordered states. One of them is the 
antiferro-quadrupolar phase (AFQ) whose on-site vector, Eq. (O, is 
mutually orthogonal to each other, and the other is the standard 120° 
AFM phase which can be characterized by calculating the inner prod- 
ucts between each pair of the different nearest-neighbor on-site vec- 
tors. 



A. Site-factorized state studies 



with # _ # running over all the bonds on the lattice; run- 
ning over all the triangles, up- and down-triangles, on the lat- 
tice, and (123) are the sites on the triangles labeled coun- 
terclockwise, Fig. 03 Pi 2 is the nearest-neighbor two-site 
exchange operator and P123 is the three-site ring exchange 
operator, which permutes the fermions on the triangles as 
Pjki \a, (3, 7) = I7, a, f3). Without the ring exchange terms, 
previous studies of the site-factorized ansatz on the triangu- 
lar lattice predicted a three-sublattice-ordered state 10-32 which 
was recently confirmed by Density Matrix Renormalization 
Group (DMRG) and infinite Projected Entangled-Pair States 
(iPEPS) analysis 2 !. 

In this section, we will focus on the SU(3) model with AFM 
three-site ring exchanges. In order to study the magnetic or- 
dered phases, we first study the model using the site-factorized 
states below in Sec. Ill Al Later in Sec. Ill Bl we will present 
our studies using the slave-fermion trial states focusing on the 
non-magnetic phases. 



In this subsection, we consider the site-factorized state* 3 - - 
defined as 

i*>=ni*i>' (2) 

with 

\Xj) = aj\x)j + b^y), + Cj\z)j, (3) 

where we fix the overall phase by setting the phase of Oj to be 
zero such that a,j E K and bj, Cj G C and |aj| 2 + \bj\ 2 + 
\cj\ 2 = 1. Above, we used the usual time-reversal invariant 
basis of the SU(3) fundamental representation^, defined as 

with \S Z = ±1) = I ± 1) and \S Z = 0) = |0). 

According to the parametrization of the \Xj) vector in 
Eq. (0 along with the constraint, at each site there are 4 in- 
dependent parameters. For a lattice with N x N sites, there 
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are AN 2 independent parameters for the site-factorized state, 
Eq. ©. We numerically find the optimized (lowest) site- 
factorized state energy, E s f = (s\Hsu(3) l s )> on a 3 x 3, and 
on a 6 x 6 triangular lattice for a certain K3 while J = 1 using 
the gradient descent method. 

The phase diagram is shown in Fig. [2] Compared with 
the site-factorized state studies on the SU(3) model with 
the FM K%, which contains the FM ordered state in addi- 
tion to the three-sublattice-ordered phase, called the antiferro- 
quadrupolar phase (AFQ) in this paper—, the phase dia- 
gram along the AFM K 3 / J-axis contains two different three- 
sublattice-ordered states, the AFQ and the standard 120° 
AFM. The AFQ is the same to the three-sublattice-ordered 
phase found in the FM K3 side,— with each on-site vector 
(\Xj e A), \XjeB), and \Xj<=c)) mutually orthogonal to each 
other. The 120° AFM can be characterized by calculating the 
inner products between each pair of the on-site vectors. 

The energies of both states can be calculated analytically. 
The site-factorized state energy of the AFQ is exactly zero ( 
E A fq = 0) and that of the AFM is E A fm = § J - \K$. 
The transition between these two states is at K3/J = 3/2, 
which is consistent with our numerical calculation. We em- 
phasize that site-factorized ansatz studies though give poor 
energies due to the ignorance of the quantum fluctuations, it 
can give us the correct order in energetics among the magnetic 
ordered phases. Comparing the result of magnetic ordered 
phases with that in Ref. |24j, we can see the order in energet- 
ics between only the magnetic ordered phases are consistent. 
For the corrections for the energies of these magnetic ordered 
phases, we can perform linear flavor wave theory calculations 
but these calculations still can not allow us to compare the 
energetics between the magnetic ordered states here and the 
non-magnetic ordered states obtained in the next subsection. 



B. Slave-fermion trial states and energetics 

In this subsection, we follow the approach similar to the one 
outlined in Ref. [33] for the S = 1 spin model. We write the 
spin operators in terms of three flavors of fermionic spinons, 



r, 



Sf 



(5) 



with a, (3, 7 € {x, y, z} and j is the site label. Rewriting 
the spin operator in terms of fermionic spinons enlarges the 
Hilbert space. To recover the physical subspace, a local con- 
straint on the fermions has to be enforced, 



1. 



(6) 



The exchange operators in terms of fermions are 

p _ f«t fP fP\ fa 

~ h h k Jk > 



(7) 



where £ Q = £ q= 



p _ f at fP fP\ fl f 7t fa 

^jkl — J j Jj Jk Jk Jl Jl ' 

a^7 



(8) 



z> V, * 



and similar relations for (3 and 7. 



The Hamiltonian, Eq. (OJ, can be re-expressed as 



H 



a 13 



a/37 



/r/f/ 2 pT / 2 7 / 3 7T /3 a +H.c. 



Below, we will calculate the trial energies using the slave 
fermion trial states. In order to check which mean-field ansatz 
states are most relevant for this model, we perform numeri- 
cally "full optimization" of the mean-field energy, Eq. ( fTTT i. on 
a triangular lattice with 100 x 100 3-site unit cells, by treating 
Xjfe-s and 8jk-$ as varying variables. In the numerical opti- 
mization, there are totally 18 variables (9 Xjk and 9 9jk) and 
we take tjk = 1, fif = /U. Numerics suggest that, at the 
mean-field level, the main competing states are the "trimer" 

state, the <J>n CSL, and the $2, CSL^i 

3 t" 

When we perform numerical calculations, we relax the con- 
straint of the fermion number for each flavor to be 



/ j-af ra\ 

\J 7" J 7 /trial 



1 



(9) 



A convenient formulation of the mean field is to consider a 
general SU(3)-rotation invariant trial Hamiltonian 



,fce 



Htnal [*J 
(jk) a 

EE/o/r/ 



(10) 



J a 



with tjk being the hopping amplitude, 6jk being the phase 
of the hopping tjk in different mean-field ansatz states, and 
fj,j being the chemical potential which can be used to satisfy 
the constraint, Eq. With the trial Hamiltonian above, we 
can find the ground state and use it as a trial wave function 
for the Hamiltonian i?s(/(3)> Eq. ([T). After performing "com- 
plete" Wick contractions and ignoring the constant pure den- 
sity terms, the trial energy can be expressed as 
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e M f =-^E E 



Xl2 



+ 



£xi2X2 3 X31 

A a 



E 



"l X23X32 + «2X3lXi3 + "MXSl + 



X13X32X21 

a/37 



H.c. 



(11) 



where we defined {Xj k T = (ff fk) trial- 

The slave-fermion trial states which conserve the transla- 
tional symmetry that we consider are the $n CSL and $2^ 
CSL, Fig. [3] There are three sublattices per unit cell for <£>» 
CSL and $ 2^ CSL and the directions of the arrows in Fig. [3] 
represent the phases associated with the fermion hoppings. 
Below we list the numerical values of the energy per site in 
the two states 



E$J n/3 = -0.8992 J - 0.8343^3, 



E l=2*/z = -0.5425 J - 0.8942^3. 



(12) 
(13) 



Besides the translationally invariant state, we also consider 
what we call the "trimer" state. Fig. |4] shows one exam- 
ple of the configuration of such a state in which the non- 
zero tjk form non-overlapping trimer covering of the lattice. 
These states break translational invariance, and any trimer 
covering produces such a state. Such states can have the 
lower Heisenberg exchange energy. The occupied bonds at- 
tain the maximal expectation value which is found analytically 
\x%\raax = n f = 1/3- Their contribution can be sufficient to 
produce the lowest total energy and such states are expected 
to be the lowest-energy states with K3 = 0. 



B«L = -J -0.5926^3. 



(14) 



Based on the energies of each trial state, Eq. (fl~2t-(fl~4l. we 
can analyze the phase diagram along K3J J-axis. At K3 = 0, 
the trimer state is the lowest energy state followed by the <E>|. 
CSL and the $2^ CSL. When K3 increases, the energy line 

of the $n CSL crosses that of the trimer state and $il CSL 

3 3 

becomes the lowest energy state at K3/J ~ 0.417. When 
K3/J keeps increasing, the energy of the $4^ will be lower 
than that of the $^ CSL when K3/J > 5.955 and becomes 
the lowest-energy state. The results are summarized in the 
mean-field phase diagram, Fig. [3] At the mean-field level, for 
each CSL, we also calculate the Chern number of the filled 
lowest band. The details will be illustrated in Sec. Ill CI In 

short, at the mean-field level, the CSL has a Chern number 

3 

C = +1 and the $ 2* CSL has a Chern number C = +2. 

3 

Before leaving this section, we want to remark that the 
trimer state is a singlet state around a triangular plaquette, and 
we can write down the exact singlet wave function in a closed 
form as 



\>P! 



a/57 



a,p,7 



(15) 



2?r w 2vr 



3 > 3 / 3 / 3. 
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(b) 



FIG. 3. Graphical illustrations of the $41 CSL and CSL. Both 
of these states have three sublattices per unit cell labelled as A, B, 
and C. The directions of the arrows represent the phases associated 
with the fermion hoppings.The arrows represent the direction of the 

gauge, (a) $n CSL. (b) CSL. 

3 3 



with a = x,y, z. With the trimer wave function, we can cal- 
culate the exact energy per site 



(16) 



In general, the exact energy of the trimer is lower than that of 
the AFQ^ and is lower than that of the AFM when K3/J < 
117/70 ~ 1.67. In order to compare the energies of the CSL 
with other states, we need to perform Gutzwiller projections 36 
on the mean-field CSL states we obtain above, which is be- 
yond the scope of this paper. 
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FIG. 4. Graphical illustration of the trimer state. In slave fermion 
picture, the fermionic spinons only hop around each triangular pla- 
quette and we can focus on each triangle separately. 
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FIG. 5. The phase diagram of the mean-field ansatzes. The ex- 
act wave function of the trimer state can be written down explicitly, 
Eq. ( 1 1 5b . and we can calculate the corresponding energy exactly. 



C. Berry curvature and Chern number 

At the mean-field level we can characterize each CSL 
discussed above by calculating its corresponding Chern 
number 37 . Full theory can have other excitations which we 
will ignore in this paper. Since each flavor of fermions have a 
filling v = 1/3 per site and forms a band insulator, there is a 
well-defined Chern number. The Chern number (which is the 
Berry phase in units of 2ir) of a many-body state at band a is 
an integral invariant in the boundary phase space, 



C a = 7T / dk 1 dk 2 
^ Jbz 



dki 



dk 2 



jig 

dk 2 



4>a 

dki 



,(17) 



where ip a is the wave-function of the band a, ki/ 2 represent 
the momentum along vectors e\ and e 2 , Fig.Q] and BZ stands 
for the first Brillouin zone. Using Eq. (|T7t , we obtain the 
Chern numbers for the lowest filled band in CSL and <f> 2jl 
CSL. We find that C = +1 for $n CSL and C = +2 for $i 

3 3 

CSL. Since the band spectra for each flavor of fermions are 
the same, we focus on one flavor of fermions. Fig. [6] shows 
the edge states of each CSL and the Chern numbers of each 
bulk band. 

It is known that the effective theory of the bulks of CSL 
are described by the effective Chern-Simons theory, which 
gives the low-energy excitations with self- and mutual- braid- 
ing statistics. Besides, the sharpest differences between dif- 
ferent CSL lie on the edges. According to the edge-bulk 
correspondence, 38 we can obtain the corresponding edge the- 
ory from the Chern-Simons description of the bulks. The ef- 
fective Chern-Simons theories hence are useful for character- 




-7T 







7T 



(a) 




-7T 



7T 



(b) 



FIG. 6. Edge states and Chern numbers of each band for a flavor of 
fermions in the trial uniform-flux states. The lowest band is com- 
pletely filled for each flavor of fermions. (a) The edge states of the 

<1?2l CSL. Each flavor of fermions show the same edge states and 

3 ° 
bulk bands whose filled lowest band is of Chern number C = +1. 

(b) The edge states of the $2^ CSL. The lowest filled bands for 

3 

three flavors of fermions show Chern numbers C = +2. The 

3 

CSL can be obtained by inserting a tt flux to each triangle along 
with flipping the arrow directions. The two spectra can be related by 

e *2,/3( k ) = - e 4V3( k + 7r )- 



izations of each CSL and will given below. 

Before giving the effective Chern-Simons theory for the 
filled lowest fermion band in the CSL, it is convenient to trans- 
form the basis from {f x J y J z } to {/+i, /-i, /o}, where 
/±i carry S z quantum number ±1 and /o carries S z quantum 
number 0. Based on Eq. ©, we can transform the fermions 
with 



V2 



(18) 



In {/±i, fa} basis, the spin operator can be represented as 

S^ = S a +iS* = v^(/| 1 /o + / '/-i) 1 (19) 



s z 



f+lf+1 - f-lf-1- 



(20) 
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In this new basis, below we will give the effective Chern- 
Simons theory for each CSL and its corresponding edge the- 
ory at the mean field. 



1. Effective Chern-Simons theory for mean-field $ w CSL 

The three flavors of fermions, {/±i, fo}, fill the lowest 
bands with Chern numbers C±\ = Co = +1. Conserved 
fermion currents J£ can be expressed in terms of dynamical 
U(l) gauge fields a™ as J£ = ^-d v a^ with m = ±1, 0, 
where summation over repeated /it, I/, A is assumed. The 
fermion band structure in this CSL is described by the 
following U 3 (l) Chern-Simons theory: 



C f = 



47T 



C m a™d v al 



711 — — 1 



The $zl CSL has the spin quantum Hall conductance (in units 

< y = = 2. (27) 

Besides, there are two different anyon excitations, both having 
the self statistical angle 9 = 2n/3. Their mutual braiding 
statistics is 6' = 2n/3. This <£n CSL is similar to the CSL 
discussed in the footnote [35] of Ref . |40[ 

Once we know the effective Chern-Simons theory for the 
mean-field CSL states, the corresponding edge theory can 
be obtained from bulk-edge correspondence^. The effective 
edge theory (assume the edge is along x-direction) for <E>il 
CSL is 

£-edge ' ^ ,,<)"■',<), 'J , - Vij d x (j)id x (j)j) 



<M E m -< 

\m— — 1 



"2= E qi(A s Q z d x <t>i 

i=l,2 



(28) 



—. — Cuaiudvdjx + — — tjA s u d v ai\, 

'ITT Z7T 



(21) 



where is the gauge potential that couple to the S z spin 
density and current, I, J = 1,2,3, and 



(22) 



The local constraint of the fermions, Eq. (|6), can be written in 
a covariant form: 














i 


1 
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• 














where K matrix and q vector are defined in Eq. (1261 1. The V 
matrix is positive-definite real symmetric which determine the 
velocities of edge modes and the number of right movers and 
left movers are determined by the number of positive and neg- 
ative eigenvalues. We note that in this case, the eigenvalues of 
K arel positive and the edges are chiral and stable, at least at 
the mean-field level here. 

The S z density on the edge are given by the defined bosons 
in the $i CSL regime as follows: 



3 1 ' 

i=l,2 



dx4>i{x) 
2ir ' 



(29) 



2?r 



rn — — 1 



2- 



d v a\, (23) 



where a M is a non-dynamical (constant) background field, 
whose density J° = (d x a y — d y a x )/2ir = R.H.S. of Eq. ((6]). 
The constraint can be implemented by introducing an extra 
U(l) gauge field as a Lagrangian multiplier, 



and (j>i carries S z quantum number +1. <p2 carries S z quan- 
tum number —1. 

We can also write down the bosonized expression of the 
transverse component of the spin on the edge: 



9+ 

°<E> - 



(x) ~ erKtyy+fa) + e »(0i+20 2 ) _ 



(30) 



constraint 



E < - ^ 

Ym=— 1 / 



~2~tV 



(E 



a/A - a\ 



(24) 



After integrating out the gauge field b^ and <Z2 M (or a° in the 
original language), we can obtain the low-energy theory of the 
$n CSL stated 



£/ 4~ ^constraint 



In 



2- 



q l Af [ Z d^a lX , (25) 



with i, j = 1,2 (where we relabel I,J = 1,3 4 i,j = 1,2) 
and the 2x2 if-matrix and q-vector are 



The edge boson fields in <E>^ CSL satisfy the Kac-Moody al- 
gebra: 

[4> j (x),d x 4> k {y)]=i27r(K- 1 ) jk 6{x-y). (31) 

At the mean-field level, because of the gapless edge excita- 
tions in the CSL phase, the transverse spin components should 
show power law (S + (x, t)S~ (0, t)) ~ | x \~ p , where pis some 
number determined by the details of the matrix V in $^ CSL. 

2. Effective Chern-Simons theory for mean-field $2. CSL 

Now the three flavors of fermions fill the lowest bands with 
Chern numbers C±\ = Co = +2. The fermion band structure 
in this $2^ CSL is described by the following U 6 (l) Chern- 
Simons theory: 

C f = ^ E E [Om+^'^Ow 



K = 













i 




-J 



(26) 



47T 

i,fc=l,2 m=±l,0 

fiivX ^yvX 



4tt 



■Cijaipdvajx 



2n 



(32) 
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where is the gauge potential that couple to the S z spin 
density and current, 7, J = 1, 2, 6, and 



C= (1)6X6, t= (1 1 -1 -1) 



(33) 



where (l) n xn represents n x n identity matrix. The local 
constraint now can be written as, 



~2n~ 



J2 d - a ^ = E J i = J " = ~^r d ^ x - (34) 



i=i 



The constraint can be implemented by introducing an extra 
U(l) gauge field b M as a Lagrangian multiplier, 



constraint 



~2n~ 



bud u f a/A - a A J . (35) 



After integrating out the gauge fields 6^ and the gauge field, 
03,^ (or a?/*)' we can obtain the effective Chern-Simons the- 
ory for the $ 2* CSL:ii 



£ 



C.S' 



£/ "1" ^constraint 



-in 



jx -t- — c^aa, (36) 

Z7T 



with i, j = 1, 2, .., 5 (where we relabel 7, J = 1, 2, 4, 5, 6 — ► 
i, j = 1,2,3,4, 5) and the TsT'-matrix and o'-vector are 



K' 



/2 1 1 1 1\ 
'12 111' 
112 11 
1112 1 
\l 1 1 1 2) 



(11 0-1 -1) . (37) 



Therefore, the $2^ CSL has the spin quantum Hall conduc- 
tance, 



q' T K'- l q' 



(38) 



There are five anyon excitations, each having the statistical 
angle 9 = 5ir/6. However, the mutual braiding statistics be- 
tween each two is 9' = 4ir/3. 

The effective edge theory for $2, CSL is 

i,j=l,..,5 

+^E«i(<^-^"^)' (39) 



where K' matrix and q' vector are defined in Eq. (|37V V' 
matrix is positive-definite real symmetric and the eigenvalues 
of K' are all positive and the edges are chiral and stable, at 
least at the mean-field level here. 

The S z density on the edge in the $2^ CSL regime is 



SS 2 , Or) 



E 4. 

i=l,2,4,5 



,d x4>'j{x) 

2tt ' 



(40) 



where (/)[, and <f>' 2 carry S z quantum number +1. <fi' 4 and <fi' 5 
carry quantum number — 1, The remaining <fi' 3 carries S z 
quantum number 0. 

The bosonized expression of the transverse component of 
the spin on the edge is 

S + (X) ~ e _ ^ 2 ^' 1+ ^2 + ^4+? i 5) 4- g'W+^+^ + ^s) Jf- 

_)_ e — '(02~*3) 4. e - l (^3-*s). (41) 

The edge boson fields satisfy the Kac-Moody algebra: 

[#(*), 0*0i(y)] - i27r(7v'- 1 ) jfe ,5(x - y). (42) 

At the mean-field level, the transverse spin components 
should show power law (S + (x, t)S~~ (0, t)) ~ | | 3=1 , where 
is some number determined by the details of the V' matrix 
in $2, CSL. 



III. DISCUSSION 

We study the SU(3) model with AFM three-site ring ex- 
changes. The site-factorized state studies over the magnetic 
ordered states show two different three-sublattice-ordered 
states, AFQ and 120° AFM, along the K 3 /J axis. We later 
use slave-fermion trial states to study the non-magnetic states. 
We find that there are three competing trial states-trimer, <E>ii 
CSL, and $2, CSL. At the mean-field level, the CSL 

has Chern number C$ „ = 1 with a finite spin quantum Hall 

3 

conductance a xy = 2. There are two different anyon excita- 
tions, both having statistical angle 9 = 2ir/3, but their mutual 
braiding statistics is 9' = 27r/3;Interestingly the $2^ CSL 
has Chern number C$ 2ll — 2 with finite spin quantum Hall 

~3~ 

conductance a xy = 4. There are five different anyon excita- 
tions, each having statistical angle 9 = 5n/6, but the mutual 
braiding statistics of each two is 9' = Att/3. 

Within the scope of this work, it is not possible to compare 
directly the energetics between the magnetic ordered states 
and the non-magnetic states since we used two different meth- 
ods to explore the phases in this model. The only wave func- 
tion of the non-magnetic state that we know is the trimer state 
whose energy can be lowered than that of the 120° AFM in 
a small parameter window as discussed in the last paragraph 
of Sec. Ill B 1 Since at the mean field, the uniform-flux CSL 
have lower energies than that of the trimer in a large parame- 
ter regime, the comparison above gives us the possibility that 
the uniform-flux CSL will be realized after performing the 
Gutzwiller projection on the mean-field states. 

It is possible to perform linear flavor wave theory (LFWT) 
as illustrated in Ref. Hj] to study the effects of quantum fluc- 
tuations on the magnetic ordered states. However, the quan- 
tum fluctuations corrections in this model actually lower the 
energies of the magnetic ordered states as shown in the note 
l35l for the AFQ. Such qualitative arguments about the quan- 
tum fluctuations destroying the magnetic order hence can not 
be directly applied here and more qualitative analysis such 
as Variational Monte Carlo (VMC), DMRG, and iPEPs are 
needed. The parameter regime in which the uniform-flux CSL 
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are more likely to be detected by the numerics can be seen as 
follows. According to the VMC results in Ref. [24j the in- 
teresting d x + idy CSL is realized in a small window around 
K3 1 J ~ 1 between the two magnetic ordered phase- AFQ 
and 120° AFM. Intuitive, in our mean field picture, since the 
uniform-flux CSL have much lower energies than those of 
the gapless QSL, which lead to the interesting d x + id y CSL 
via plausible pairing mechanisms, the small window of the 
d x + idy CSL is likely replaced by that of the uniform-flux 
CSL, especially the $2: CSL since $2^ CSL is deep inside 
the 120° AFM regime.' 

In the present work, we ignore the gauge fluctuations in 
the mean-field slave-fermion analysis. The gauge fluctuations 
around our mean-field results can be described by the com- 
pact U(l) lattice gauge theory in (2+l)D, which brings out 
the confinement issue due to the proliferation of the topo- 
logical defects. Indeed, in (2+l)D, if the partons (spinons) 
in the U(l) slave-fermion theory are gapped without break- 
ing Time-Reversal symmetry (such as 7r-flux) and the filled 
bands are trivial, a.k.a. Chern number C = 0, the topolog- 
ical defects along the space-time known as the instantons or 
Dirac monopoles are proliferating to induce the instabilities 
that lead the gapped QSL to some order states which break 
translational invariance.— In sharp contrast, if the filled par- 
ton bands are topological as in the present uniform-flux CSL, 
the effective theory of the gauge fluctuations is described by 
the U(l) lattice gauge theory with a Chern-Simons term. The 
Chern-Simons term has dramatic effects on the behavior of 
the Dirac monopoles. The Dirac monopoles become charged 
and linearly confined due to the presence of a Chern-Simons 
ter m 3-4 , and the uniform-flux CSL can remain stable against 
gauge fluctuations. 

The topological order in the CSL phases can be par- 
tially detected using entanglement spectrum and the topo- 
logical entanglement entropy (Renyi entropy or von Neu- 
mann entropy)4£^ The Chem number of the CSL phases 
can be possibly extracted using Exact Diagonalization^ Fur- 
thermore, recently it is also suggested that the quasi-particle 



self and mutual braiding statistics can be extracted from the 
ground states with minimum entropies using VMC^i 

For a possible future direction of the present theoretical 
model, even though it may be difficult to detect the inter- 
esting uniform-flux CSL in the cold atom experiments, since 
closer to the metal-Mott insulator phase transition the higher 
order terms such as four-site ring exchange terms also needs 
consideration, the present model can serve as a theoretical 
starting point for a lattice model realization of the interest- 
ing E/s*(l) symmetry protected topological phase (SPT) in 
the featureless spin- 1 frustrated magnet4£ Following the logic 
detailed in Ref. |4CJ the Us* (1) SPT can be possibly realized 
in the present model by introducing a Time-Reveral breaking 
(TRB) and S z -conserving two-site exchange term. The possi- 
ble Hamiltonian is 

H u{1 )spt = jJ2 Pl2 + K * Yl t Pl23 + H - c -l + 

~ A 

m-m^P 

where a, f3 = x, y, and z. The introduced TRB and S z - 
conserving P^f in the slave-fermion representation is ex- 
pressed as Pf/ = eWv-ofof^fff^fg, which contains 

the complex phases as e l ^ e ° 2 ~ e ^ with 9f 2 = 6f 2 = —#12 
that realizes the Us- (I) symmetry of the rotation about S z - 
axis. This Ktrb term adds another axis to the parameter 
spaces in the present model and whether or not the U~s~(l) 
SPT can be energetically favored is interesting by its own. 
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